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1 . 0  Introduction 


The  Atlantic  Hurricane  Wind  Probability  (WINDPA)  model  is  a 
continuation  of  the  Atlantic  Hurricane  Strike  Probability  Pro¬ 
gram  developed  by  Jarrell*  (1981).  This  model  extends  the 
Atlantic  strike  probability  model  using  wind  probability  con¬ 
cepts  similar  to  those  developed  and  presented  for  the  western 
Pacific  Ocean  basin  (Jarrell  ,  1981).  Those  concepts  will  not 
be  presented  again  here  except  to  illustrate  the  differences  in 
methodology  between  the  two  models.  The  parameterization  of  the 
asymmetric  wind  distribution  around  the  tropical  cyclone  and  the 
inference  of  30  and  50  knot  wind  radii  are  different  and  will  be 
described. 

2.0  Model  Description 

2.1  Atlantic  Hurricane  Strike  Probability  Program 

The  strike  probability  concept  was  developed  by  Jarrell3 
(1978)  and  has  been  successfully  applied  to  the  eastern  and 


1 Jarrell,  J.D.,1981:  Atlantic  Hurricane  Strike  Probability 
Program;  NAVENVPREDRSCHFAC  Contractor  Report  CR81-04. 

2 Jarrell,  J.D.,1981:  Tropical  Cyclone  Wind  Probability 
Forecasting  (WINDP) ; NAVENVPREDRSCHFAC  Contractor  Report  CR81-03. 

3Jarrell,  J.D.,1978:  Tropical  Cyclone  Strike  Probability 
Forecasting;  NAVENVPREDRSCHFAC  Contractor  Report  CR78-01. 
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western  Pacific  as  well  as  the  Atlantic  basin.  The  development 
of  strike  probability  theory  and  its  extension  tc  wind 
probability  is  based  upon  three  basic  assumptions: 

(a)  All  forecasts  related  to  tropical  cyclones  are 
subject  to  some  degree  of  error. 

(b)  The  size  of  the  position  forecast  error  is 
statistically  related  to  the  relative  difficulty 
of  the  forecast. 

(c)  The  occurrence  of  position  forecast  errors  is 
random  and  approximates  a  bivariate  normal  prob¬ 
ability  distribution  ;  normal  in  both  N-s  and  E-w 
directions . 

Investigators  verified  these  assumptions  in  each  basin  in 
separate  studies;  Nicklin* ( 1977 )  in  the  western  Pacific, 
Thompson  and  Elsberry5( 1979 )  in  the  eastern  Pacific  and 
Crutcher6( 1980)  in  the  Atlantic.  In  these  studies  statistical 


4Nicklin,  D.S.,  1977:  A  Statistical  Analysis  of  Western 
Pacific  Tropical  Cyclone  Forecast  Errors;  Naval  Postgraduate 
School;  M.S.  Thesis. 

5Thcmpson,  W.J.  and  R.L.  Elsberry,  1979:  A  Statistical 
Analysis  of  Eastern  Pacific  Tropical  Cyclone  Forecast  Errors; 
12th  Tech.  Conf.  on  Hurricanes  and  Tropical  Meteorology;  New 
Orleans;  April,  1979. 

6Crutcher,  H.L.,  1980:  Tropical  Storm  Forecast  Error  and 
the  Bivariate  Normal  Distribution.  13th  Tech.  Conf.  on 
Hurricanes  and  Tropical  Meteorology;  AMS;  Miami,  FL;  Dec, 
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methods  were  used  to  group  the  forecasts  into  three  classes  of 
relative  forecast  difficulty  of  easy,  average  and  difficult. 

The  resultant  classifications  yielded  respective  relative  fore¬ 
cast  error  groups  of  below  average  errors  (Class  I),  average 
errors  (Class  II)  and  larger  than  average  errors  (Class  lit). 
Statistical  parameters  were  developed  in  each  study  to  describe 
the  bivariate  normal  distributions  for  each  of  these  three 
classes  and  for  24,  48  and  72  hour  forecasts.  Although  relative 
separation  of  groups  varied  in  the  three  studies,  three  distinct 
groups  emerged  in  each. 

The  statistical  data  provided  the  basis  for  a  strike  orob- 
ability  model  for  the  respective  basins.  Using  forecast  posi¬ 
tions  and  integrating  in  time  and  space,  probabilities  of  a 
tropical  cyclone  strike  are  derived  for  forecasts  out  to  72 
hours. 

Crutcher  (1980)  used  a  clustering  model  (NORMIX)  to  develop 
three  discrete  bivariate  normal  populations  of  Atlantic  position 
forecast  errors.  Because  uncertainty  is  involved  in  catego¬ 
rizing  a  particular  forecast  as  a  member  of  a  single  forecast 
difficulty  class  (even  after-the-fact  it  is  a  probabilistic 
problem),  a  method  (Jarrell1,  1981)  was  developed  to  assign 
"probability  of  class  membership"  to  each  forecast.  A  series  of 
stepwise  linear  regression  equations  were  fit  to  predictors 
relating  class  probabilities  to  the  predictors.  A  forecast 
would  have  a  likelihood  of  belonging  to  Class  I,  II  or  III  and 
possessing  those  populations  characteristics.  Each  of 


these  predicted  probabilities  is  constrained  to  0  <  Px  <  1 
and  that  Pj  +  P2  +  P3  =  1.  The  actual  strike  probability 
then  is  derived  from  three  separate  runs;  each  run  using  the 
bivariate  parameters  from  a  different  population.  The  final 
probabilities  are  a  sum  weighted  by  the  predicted  class 
probabilities.  This  method  also  forms  the  basis  of  the  Atlantic 
wind  probability  model  class  selection. 

2.2  Western  Pacific  Wind  Probability  Model  (WINDP) 

The  wind  probability  model  enhances  the  strike  model  by 
providing  probabilities  of  sustained  winds  of  30  and  50  knots  at 
a  user  specified  site.  A  definite  advantage  is  realized  for  the 
user  in  that  he  can  relate  probabilities  to  a  recognized 
destructive  wind  force  (i.e.,  30  or  50  knots)  rather  than  to  an 
arbitrarily  selected  distance  which  constitutes  a  tropical 
cyclone  strike.  The  WINDP  model  is  based  on  the  strike  model 
plus  an  evaluation  of  maximum  wind  and  wind  radius  forecasts. 
This  evaluation  is  based  on  the  Riehl7  (1963)  profile 
VR1/2  =  constant,  or  Rc  =  (Vm/Vc)2  Rm'  where 
Rc,  Vc  are  a  radius  and  wind  speed  of  concern  and  Vm,  Rm 
are  a  maximum  wind  speed  and  associated  radius.  Using  forecast 
warning  data  to  obtain  Vm  and  inferring  Rm  from  the  radius 
of  50  or  30  knot  winds,  Rc  can  then  be  approximated. 

7Riehl,  H. ,  1963:  Some  Relations  Between  Wind  and  Thermal 
Structure  of  Steady  State  Hurricanes;  Journal  of  Atmospheric 
Science;  Vol.  20;  July,  1963,  pp  276-287. 
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Rm  and  Vm  are  treated  as  random  variables.  They  are 
forecast  values  and  are  assumed  to  be  related.  All  values  of 
Vm  are  considered  and  its  probability  of  occurrence  estimated 
from  a  wind  error  algorithm.  For  each  such  Vm,  an  expected 
value  of  Rm  is  predicted  and  from  Vr^2  =  constant,  a 
value  of  Rc  is  computed.  The  marginal  probability  of  the 
point  of  interest  receiving  winds  >  Vc  (given  the  currently 
assigned  Vm  value)  is  the  probability  that  the  cyclone  passes 
within  distance  Rc  of  the  point.  The  total  probability  that  a 
point  will  receive  winds  of  at  least  Vc  at  a  time  step  is  then 
estimated  by  summing  over  all  Vm  values. 

Asymmetry  of  tropical  cyclone  storm  pattern  (i.e.,  a  larger 
semicircle  for  equal  isotach  winds  to  the  right  of  storm  track) 
is  treated  as  a  series  of  circular  isotachs  offset  to  the  right 
of  the  forecast  track.  The  forward  speed  of  motion  (S)  of  a 
storm  is  added  to  a  stationary  storm  maximum  wind  (Vg)  to 
approximate  the  maximum  wind  in  right  semicircle (Vm=Vs+S ) , 
and  subtracted  (Vs-S)to  approximate  the  left  semicircle.  An 
offset  distance  is  then  calculated  by  substituting  into  the 
Riehl  equation  and  simplifying  to  obtain: 

D  =  (2S(Vm-S)/Vc2)Rm  , 

where  V  =  V_+S.  The  offset  is  then  applied  to  the  left  of 
m  s 

the  point  of  concern  (to  compensate  for  larger  storm  semicircle 
on  right  side  of  storm)  and  used  to  determine  a  center  for 
integration.  Both  wind  profile  derivation  and  asymmetric  offset 
are  fully  developed  by  Jarrell2 ( 1981 )  for  the  interested 
reader . 
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2.3  Atlantic  Wind  Probability  Model 

The  basic  differences  to  be  observed  in  the  Atlantic  model 
as  compared  to  the  western  Pacific  are  minor  in  nature  and  lie 
in  the  derivation  of  the  50  and  30  knot  wind  profiles  and  the 
method  in  which  asymmetry  is  handled. 

The  basis  for  derivation  of  the  wind  profiles  is  from  work 
by  Tsui8(1980).  Tsui  extracted  wind  radius  data  from  tropical 
cyclone  warnings  issued  by  the  Joint  Typhoon  Warning  Center  on 
Guam  for  a  12-year  period  (1966  to  1977).  From  this  data  set  he 
estimated  the  profile  of  the  tangential  wind  speed  along  the 
radial  axis  to  be  exponential.  The  study  also  indicated  that 
the  size  of  a  tropical  cyclone  is  statistically  related  to 
maximum  wind  and  persistence  and  that  the  asymmetric  shape  of 
the  isotach  is  correlated  to  the  speed  of  movement. 

The  asymmetric  shape  of  the  wind  distribution  was 
parameterized  by  expressing  all  directional  radii  as  fractions 
of  the  radius  on  the  right  side.  For  example  the  average 
fraction  of  the  left  side  of  the  50-knot  wind  radius  is  0.73  and 
the  30-knot  wind  radius  is  0.81.  This  asymmetry  is  related  to 
speed  of  storm  movement  and  is  more  pronounced  at  higher 
speeds. 


®Tsui,  T.L.,  1980:  Surface  Wind  Distribution  of  Western 
North  Pacific  Tropical  Cylones;  13th  Tech.  Conf.  on  Hurricanes 
and  Tropical  Meteorology;  Miami,  FL;  Dec,  1980. 
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By  scaling  of  both  the  wind  and  the  right  side  wind  radius, 
data  was  then  composited  and  fitted  to  a  single  profile, 

v/vmax  =  exp(-0.693R>  , 

where  V  is  the  wind  speed  of  interest,  Vmax  is  the  maximum 
wind  (intensity)  and  R  is  the  ratio  of  the  radius  associated 
with  V  to  the  radius  associated  with  winds  of  one  half 
vmax  *n  a  tropical  cyclone. 

The  advantage  of  the  Tsui  profile  over  the  Riehl  profile  for 
this  purpose  relates  to  the  scaling  radius  used.  Riehl  uses  the 
radius  of  maximum  winds  (Rm)  as  his  scaling  radius.  Rm  is 
small  (on  the  order  of  20-30  miles)  and  is  subject  to  very  large 
percentage  errors.  Tsui,  on  the  other  hand,  uses  the  radius  of 
50%  of  the  maximum  wind  (R^alf*  as  his  scaling  radius. 

Rhalf  lar<3er  (typically  50-150  miles)  and  is  also  far 
enough  removed  from  the  central  core  of  winds  to  be  identifiable 
in  synoptic  ship  reports.  For  both  of  these  reasons, 

Rhaif  is  subject  to  smaller  percentaqe  errors  than  Rm. 

Since  in  both  cases  the  scaling  radius  is  used  in  ratio  to  some 
other  isotach  radius,  percentage  error  is  the  relevant  measure 
of  accuracy. 


With  the  above  relationship  plus  knowledge  of  the  maximum 
wind  and  its  radius  and  one  observed  wind,  an  estimate  of  the 
wind  field  around  a  tropical  cyclone  can  be  deduced.  Thus  with 
Tsui's  wind  profile  any  wind  radii  may  be  estimated  and  storm 
asymmetry  handled  with  a  simple  empirical  relationship. 

3.0  Testing  the  Atlantic  Wind  Probability  Program  (WINDPA) 

The  methodology  used  in  testing  WINDPA  predicted  values 
against  observed  values  is  identical  to  that  used  in  the  western 
Pacific  by  Jarrell2 ( 1981 ) .  An  array  of  30  points  in  the 
Atlantic  and  Gulf  of  Mexico  was  selected  (figure  1).  WINDPA 
values  for  30  and  50  knots  were  calculated  at  12  hour  intervals 
from  the  effective  synoptic  time  of  the  National  Hurricane 
Center  (NHC)  forecasts  for  the  1980  season.  Since  most  of  these 
30  points  are  not  observing  stations,  actual  verifying  winds 
were  not  generally  available.  Consequently  a  verifying  "warning 
time"  probability  greater  than  50%  constituted  a  verifying 
strike . 

Tables  1,  2,  3  and  4  compare  the  expected  to  the  observed 
occurrences  of  30  and  50  knot  winds.  Predictions  are 
associated  with  percentage  groups  of  increasing  width,  <  1/2%, 
1/2  to  1  1/2%,  .  .  .  etc.  Time  integrated  probabilities  were 
verified  only  if  a  continuous  record  was  available  over  the 
entire  time  period. 
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80 

9 

5 

154-314 

25 

6 

5 

25 

5 

4 

25 

5 

7 

314-634 

16 

7 

4' 

12 

5 

5 

6 

2 

2 

>634 

7 

5 

5 

4 

3 

3 

3 

3 

3 

ALL 

Table  5- 

5160 

Time 

23  15  3000  24 

Integrated  STRIKPA. 

16 

1860 

28 

17* 

INST 

24 

Hr 

48 

Hr 

72 

Hr 

A<P<B 

N 

E 

0 

N 

E 

O 

N 

E 

O 

<  4% 

5707 

0 

0 

3684 

0 

0 

2727 

1 

1 

4-  14 

35 

0 

0 

120 

1 

0 

242 

2 

2 

14-  34 

22 

1 

0 

77 

2 

3 

71 

1 

2 

34-  74 

21 

1 

0 

19 

1 

1 

0 

0 

0 

74-154 

17 

2 

2 

0 

0 

0 

0 

0 

0 

154-314 

18 

4 

3 

0 

0 

0 

0 

0 

0 

314-634 

0 

0 

0 

0 

0 

0 

0 

0 

0 

>634 

0 

0 

0 

0 

0 

0 

0 

0 

0 

ALL 

Table  6. 

5820  8  5  3900 

Instantaneous  STRIKPA. 

4 

4 

3030 

4 

5 

*difference  significant  at  5%  level.  (None  of  these  are  significant 
when  it  is  assumed  that  only  1/3  of  cases  are  independent.) 
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Significance  of  the  differences  between  the  expected  and  the 
observed,  as  discussed  in  previous  reports  ( Jarrell2, 1981) ,  is 
difficult  to  assess,  but  using  previously  developed  tests, 
agreement  appears  to  be  very  good.  The  instantaneous  prob¬ 
abilities  show  excellent  correlation  at  all  forecast  lengths 
with  minor  underforecasting  at  the  longer  forecast  times  (48  and 
72  hr).  The  time  integrated  probabilities  displayed  a  slight 
tendency  toward  overforecasting  for  both  50  kt  and  30  kt  winds. 
This  was  not  at  a  significant  level  (see  note  accompanying 
tables)  and  the  overall  results  are  considered  statistically 
sound.  An  instantaneous  and  a  time  integrated  strike 
probability  (STRIKPA)  was  run  on  the  same  data  base  (Table  5  and 
6)  and  showed  comparable  agreement  between  predicted  and 
observed  occurrences. 

4.0  Operational  Products 

The  Atlantic  wind  probability  program  will  be  available  for 
the  preselected  points  receiving  the  Atlantic  strike  program. 

(The  STRIKPA  program  is  actually  integrated  into  the  WINDPA  pro¬ 
gram).  Probabilities  will  be  given  in  two  modes,  instantaneous 
and  time  integrated,  and  at  0,  12,  24,  36,  48,  60  and  72  hours 
after  the  warning  time.  The  instantaneous  probability  will  be 
the  probability  at  the  stated  time  (i.e.,  12  hr)  and  the  time 
integrated  will  be  summed  for  the  0  to  X  hour  time  interval  for 
an  estimate  of  the  probability  that  the  event  will  be  observed 
within  that  period  cf  time. 
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The  crreatest  source  of  probable  error  for  the  WINDPA  program 


will  be  erroneous  input  data.  An  internal  check  for  unusual 
motion  (expected  to  occur  only  5%  of  the  time  in  nature)  will  be 
made  and  suspect  motion  flagged.  The  user  should  then  recheck 
input  data  for  accuracy. 

When  the  forecast  track  approaches  land  mass  the  forecaster 
should  be  aware  of  program  bias.  This  should  be  minor  for  sea¬ 
ward  approach  to  low  coastal  areas  or  over  smaller  islands. 
However,  in  other  cases  land  influences  will  appear  as  rapid  de¬ 
creases  in  the  instantaneous  wind  probabilities  (especially  50 
kt  winds)  near  forecast  landfall  time.  This  will  bias  prob¬ 
abilities  -  overstate  them  for  inland  sites  and  understate  them 
for  coastal  sites.  Time  integrated  probabilities  will  be  less 
biased.  This  problem  is  caused  by  wind  forecasts  being  in¬ 
fluenced  by  track  forecasts  where  landfall  is  concerned.  A  bad 
track  forecast  may  cause  a  bad  wind  forecast.  This  was  not 
accounted  for  either  in  development  nor  testing;  hence  the  test 
results  simulate  expected  actual  operational  results  and  some  of 
the  minor  disparities  between  expected  and  observed  occurrences 
no  doubt  stems  from  this  problem. 
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SUMMARY 


The  wind  probability  model  for  the  Atlantic  is  largely  based 
on  the  strike  probability  program  for  the  Atlantic  and  the  wind 
probability  program  for  the  western  Pacific  -  both  currently 
operational.  The  single  important  improvement  represented  in 
the  concept  is  the  addition  of  the  Tsui  wind  profile.  Test 
results  of  the  Atlantic  WINDPA  program  demonstrated  excellent 
agreement  between  expected  and  observed  results. 
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NATIONAL  OCEANIC  A  ATMOS.  ADHIN. 
OCEANOGRAPHIC  SERV.  DIV. 

6010  EXECUTIVE  BLVO. 

ROCKVILLE,  HO  20852 

DIRECTOR,  PROGRAMS  OFFICE,  RX3 
NOAA  RSCH.  LAB. 

BOULOER,  CO  80302 

DIRECTOR 

NATIONAL  HURRICANE  CENTER,  NOAA 
GABLES  ONE  TOWER 
1320  S.  DIXIE  HWY. 

CORAL  GABLES,  FL  33146 

NATIONAL  WEATHER  SERVICE 
WORLD  WEATHER  BLDG. 

ROOM  307 
5200  AUTH  ROAD 
CAMP  SPRINGS,  MO  20023 

NATIONAL  CLIMATIC  CENTER 
ATTN:  L.  PRESTON  0542X2 
FEDERAL  BLOG.  -  LIBRARY 
ASHEVILLE,  NC  2G80I 

DIRECTOR 

NATIONAL  OCEANO.  DATA  CENTER 
DrPT.  OF  COHMERi. F  .  NOAA 
ROCKVILLE,  HO  2U..52 

NATIONAL  WEATHER  SERVICE 
A1TN:  WFU.  f As! CRN  REGION 
6H6  STEWART  AVE. 

CARDEN  CUT.  Hi  115.10 


CH3EF ,  SCIENTIFIC  SERVICES 
NWS  ,  SOUTliri'H  REGION 
NOAA ,  ROOM  1 OE09 
Bit  TAYLOR  ST. 
n,  WORT  FI .  TX  76102 

MU80R010GIST  IN  CHARGE 
NWS  FORECAST  OFFICE 
UtHNCEOGY  II,  NOAA 
Ml V  YORK  UNIVERSITY 
BRONX,  NY  10453 

D2KEC1  OR 

ATLANTIC  OCEANO.  6  METEOR.  LABS. 

15  RICKENBACKER  CAUSEWAY 
VIRGINIA  KEY 
MIAMI,  EL  33149 

DIRECTOR 

ATLANTIC  MARINE  CENTER 
COAST  6  GEODETIC  SURVEY,  NOAA 
4»  WEST  YORK  STREET 
NORFOLK ,  VA  23510 

ST-NOAA  LIAISON  MANAGER 
NASA-JOHNSON  SPACE  CENTER 
ATTN :  DR.  M.  HELFERT 
HOUSTON,  TX  77058 

DIRECTOR 

INTERNATIONAL  AFFAIRS  OFFICE,  NOAA 
ATTN :  MR.  N.  JOHNSON 
6070  EXECUTIVE  BLVO. 

R0UV1LLE,  MD  20852 

DIRECTOR,  NSF 

Dll.  OF  ATMOSPHERIC  SCIENCES 
ROOM  644 

16S0  G  STREET,  NW 
WASHINGTON,  DC  20550 

LAJORATORY  FOR  ATMOS.  SCIENCES 
NASA  GOODARO  SPACE  FLIGHT  CENTER 
GNEENBFLT,  MD  20771 

PRELIMINARY  SYS.  DESIGN  GROUP 
NASA  GODDARD  SPACE  FLIGHT  CENTER 
GMENBELT,  MO  20771 

EMCUTIVE  SECRETARY.  CAO 

SUBCOMMITTEE  ON  ATMOS.  SCIENCES 
NATIONAL  SCIENCE  FOUNDATION 
ROOM  510 

1*®0  G  STREET,  NH 
WASHINGTON,  OC  2OSS0 

FEAERAL  EMERGENCY  MANAGEMENT 
AGENCY  (FIMA) 

17:5  1  STREET,  NW 
WASHINGTON.  OC  20472 

COLORADO  STATE  UNIVERSITY 
ATMOS .  SCIENCES  DEPT. 

ATTN:  LIBRARIAN 
FI.  COLLINS,  CO  80521 

CHAIRMAN 

Pf Nil  STATE  UNIVERSITY 
METEOROLOGY  DEPT. 

503  Of  IKE  BLOG. 

UNIVERSITY  PARK,  PA  16802 

CHAIRMAN 

MASSACHUSETTS  INST.  OF  TECH. 
METEOROLOGY  DEPT. 

CAMBRIDGE,  MA  02139 

DUFCTCR 

UCLA,  INST1.  OF  GEOPHYSICS 
LOS  ANGELES .  CA  90024 

■jMvn.'SITY  OF  HASHINGTON 
ATMOS  .  SC  1 1  NCI  '..  I>i  PT  . 

SI  A  1 1 1  (  ,  VIA  H..r<5 


FLORIDA  STATE  UNIVERSITY 
ENVIRONMENTAL  SCIENCES  OEPT. 
TALLAHASSEE,  FL  32306 

UNIVERSITY  OF  HAWAII 
METEOROLOGY  OEPT. 

2525  CORREA  ROAD 
HONOLULU,  HI  96822 

CHAIRMAN 

UNIVERSITY  OF  WISCONSIN 
METEOROLOGY  OEPT. 

1225  W.  DAYTON  STREET 
MADISON,  WI  53706 

TEXAS  ASM  UNIVERSITY 
METEOROLOGY  OEPT 
COLLEGE  STATION,  TX  7863 

CHAIRMAN 

RUTGERS  UNIVERSITY 
OEPT.  OF  METEOR.  S  PHYS.  OCEANO. 
COOK  COLLEGE,  P.O.  BOX  231 
NEW  BRUNSWICK,  NJ  08903 

DIRECTOR  OF  RESEARCH 
ST.  THOMAS  UNIVERSITY 
INSTITUTE  FOR  STORM  RESEARCH 
38)2  MONTROSE  BLVD. 

HOUSTON,  TX  77006 

CHAIRMAN 

SAN  JOSE  STATE  UNIVERSITY 
METEOROLOGY  DEPT. 

SAN  JOSE.  CA  95192 

SCRIPPS  INSTITUTION  OF  OCEANO. 
LIBRARY  DOCUMENTS/REPORTS  SECTION 
LA  JOLLA.  CA  92037 

DIRECTOR 

OLD  DOMINION  UNIVERSITY 
OCEANOGRAPHIC  INSTITUTE 
NORFOLK,  VA  23508 

OFFICER  IN  CHARGE 
NAVOCEAHCOMDEI 
U.S.  NAVAL  AIR  STATION 
FPO  NEW  YORK  09560 


UNIVERSITY  OF  MIAMI 
RSMAS  LIBRARY 

4600  RICKENUACKER  CAUSEWAY 
VIRCINIA  KEY 
MIAMI,  FL  33)49 

DIRECTOR 

LOUISIANA  STATE  UNIVERSITY 
ATTN:  0.  HUH,  COASTAL  STUDIES 
BATON  ROUGE,  LA  70803 

CHAIRMAN 

ATMOS.  SCIENCES  DEPT. 
UNIVERSITY  OF  VIRGINIA 
CHARLOTTESVILLE,  VA  22903 

UCLA 

ATMOSPHERIC  SCIENCES  DEPT. 

405  HILGARD  AVE. 

LOS  ANGELES,  CA  90024 

COLORADO  STATE  UNIVERSITY 
ATMOS.  SCI.  OEPT.  LIBRARY 
FOOTHILLS  CAMPUS 
FT.  COLLINS,  CO  80523 

UNIVERSITY  OF  MV  RYE AND 
METEOROLOGY  OEPT. 

COLLEGE  PARK,  MD  20742 

THE  EXECUTIVE  DIRECTOR 
AMERICAN  METEOR.  SOCIETY 
45  BEACON  STREET 
BOSTON,  HA  021  OB 

AMERICAN  MET.  SOCIETY 
METEOR.  X  GEOASTRO.  ABSTRACTS 
P.O.  BOX  1736 
WASHINGTON,  DC  20013 

DIRECTOR,  JTWC 
BOX  17 

FPO  SAN  FRANCISCO  96630 

WORLD  METEORO.  ORGANIZATDN 
ATS  01V.  (ATTN:  N .  SUZUKI) 
CH-1211,  GENEVA  20 
SWITZERLAND 


LIBRARIAN 

UNIVERSITY  01  MELBOURNE 
METEOROLOGY  OH-T. 

PARKVILll,  VICTORIA  3052 
AUSTRALIA 

BUREAU  OF  METEOROLOGY 
ATTN;  LIBRARY 
BOX  1289k,  GPO 
MELBOURNE,  VIC,  3001 
AUSTRALIA 

L I BRARY 

ATMOSPHERIC  ENVIRONMENT  SERVICE 
490S  0UFFCR1N  STREET 
DOWNSVIEW  M3II  ST 4 
ONTARIO,  CANADA 

DIRECTOR  OF  NAVAL  OCEANOGRAPHY 
A  METEOROLOGY 
MINISTRY  OF  DEFENCE 
OLD  WAR  OFFICE  BLDG. 

LONDON,  S.W.l .  ENGLAND 

METEOROLOGICAL  OFFICE  LIBRARY 

LONDON  ROAD 

BRACKNELL,  BERKSHIRE 

RG  12  2SZ 

ENGLAND 

COMMANDER  IN  CHIEF  FLEET 
ATTN:  STAFF  METEORO.  i 
OCEANOGRAPHY  OFFICER 
HORTHWOOD 
MIOOLESEX  HAG  3HP 
ENGLANO 

THE  DIRECTOR 

INDIAN  INSTITUTE  OF  TROP.  METEO. 
RAMDURG  HOUSE 
PUNE  411-005 
INDIA 

FEDERAL  COORDINATOR  FOR  MEUORO. 

SERVICES  «  SUPPORTING  RESEARCH 
6010  EXECUTIVE  6LV0. 

ROCKVILLE,  MD  20852 
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